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£S| Abstract 

o ■ 

£N| ■ We calculate the emergent spectrum of microscopic black holes, which emit copious 

amounts of thermal "Hawking" radiation, taking into account the proposition that 
(contrary to previous models) emitted quarks and gluons do not directly fragment into 
hadrons, but rather interact and form a photosphere and decrease in energy before 
fragmenting. The resulting spectrum emits copious amount of photons at energies 
JL | around lOOMeV. We find that the limit on the average universal density of black holes 

is not significantly affected by the photosphere. However we also find that gamma ray 
satellites such as EGRET and GLAST are well suited to look for nearby black holes out 
to a distance on the order of 0.3 parsecs, and conclude that if black holes are clustered 



locally as much as luminous matter, they may be directly detectable. 



submitted to Physical Review Letters 



I Introduction 

Since Hawking [||, [| first proposed that a black hole emits thermal radiation with an emission 
rate inversely proportional to its mass, there have been several calculations of the emergent 
spectrum in order to verify or at very least constrain the presence of the smallest, hence most 
luminous ones. Nominally, aside from particle spin effects and gravitational backscattering 
effects close to the black hole, calculated by Page [Q], one would expect the the emergent 
spectrum to be thermal, since Hawking showed that the black hole can be thought of as a 
black body with a temperature T = (SttGM)" 1 , where M is the mass of the black hole, G 
is Newton's gravitational constant, and we set the Boltzmann constant, k = 1 (along with 
c = 1 and h = 1). 

However, when a detailed model of the physics of the emitted particles in considered, the 
emergent spectrum becomes more complicated. For example, MacGibbon and Weber Q and 
Halzen et al. ||] have considered black holes with temperatures greater than the characteristic 
QCD energy scale A ~ 200MeV, where the black hole begins to emit quarks and gluons. 
They propose that the emitted quarks and gluons fragment into hadrons, which further 
decay into photons, electrons, neutrinos etc., and they convolve a jet code with the Hawking 
thermal spectrum to determine the emergent spectrum. For black hole temperatures above 
A, the QCD degrees of freedom dominate in the standard model and the ultimate products of 
quark and gluon fragmentation will dominate the spectrum, thus understanding the physics 
of the quarks and gluons is important for determining the spectrum. 

In this letter, we reconsider an important initial assumption made by MacGibbon and 
Weber, namely that the quarks and gluons emitted from the black hole directly fragment into 
hadrons. Rather, since the density of emitted particles around the black hole can be very 
high (much higher than nuclear density) for T > A, we propose that the quarks and gluons 
propagate through this dense plasma and lose energy via QCD bremsstrahlung and pair 
production interactions until the density of the outward-propagating plasma becomes low 
enough that the quarks finally fragment into hadrons. One can then calculate an emergent 
spectrum by convolving this collection of lower energy quarks with jet codes or fragmentation 
functions, but this spectrum will be very different than the one obtained using fragmentation 
of quarks coming directly from the black hole. In addition to determining the emergent 
spectrum of the black hole including the effect of interactions, we will also consider several 
observational consequences of the this emergent spectrum, using the EGRET and GLAST 
satellites as exemplary detectors. 

Since it is reasonable to assume that fundamental modes will be present in a thermal 
bath, let us assume that once above a temperature T > A, a black hole emits individual 
quarks and gluons. Hawking showed that the emission rate spectrum for each particle degree 
of freedom is dN/dE = T s /2-n'(ex i p(E/T) — (— l) s ) _1 , for particles of energy E, where T s is 
the absorption coefficient which in general depends upon the spin of the particle s, M, and E 
Hi, and we have assumed the black hole to be uncharged and non-rotating. In the relativistic 
limit T ^> m the total rate at which particles are directly emitted from the black hole can 
be expressed in the form A^ tot = (10~ 2 77)T, where rj is of order of the number of emitted 
relativistic particles at temperature T and can be calculated numerically [|, |J. Since this 
is the flux of particles crossing the Schwartzschild radius r s = (47rT) _1 , the density n(r) of 
emitted particles at a radius r from the black hole is then n = N /(Aitr 2 ). Expressing this in 



a more illuminating form, we obtain 

For the QCD, 77 ~ 20, thus the average particle separation, defined as d = n -1 / 3 is then 
d{r)-T-\r s /rfl\ 

QCD is an asymptotically free theory of interactions between quarks and gluons. In gen- 
eral, when particles scatter, the momentum exchanged must be at least of order their inverse 
separation d~ l . As a consequence, at distances smaller than A -1 the QCD interaction is per- 
turbative, while at larger distances, the coupling becomes so large that vacuum polarization 
and fragmentation is dominant. In particular, vacuum fragmentation of quarks and gluons 
will occur when they are separated by a distance greater than A -1 . We see from the above 
formulas that once T > A, the quarks and gluons around the black hole will not immediately 
vacuum fragment into hadrons because they are closely spaced in a kind of plasma. Rather, 
they will propagate for some distance in the dense quark-gluon plasma until the plasma 
becomes tenuous enough that vacuum fragmentation will occur. This is important because 
as the quarks and gluons propagate away from the black hole in the dense plasma, they will 
interact with each other via bremsstrahlung and pair production and decrease in energy. 

This can be seen by following the arguments of Heckler ||, who has shown that QED and 
QCD bremsstrahlung and pair production interactions become important among particles 
emitted from a black hole above some critical temperature. The essential argument stems 
from the fact that the relativistic bremsstrahlung cross section is independent of energy, and 
since the density of particles around the black hole increases with temperature (eq. [I]), there 
is a critical temperature for which the optical depth of an emitted particle becomes unity. 
At this point particles begin to scatter copiously, and a kind of photosphere forms around 
the black hole. The photosphere is a kind of fireball in the sense that the nearly thermalized 
plasma propagates outward, decreasing in temperature until i) the electrons and positrons 
annihilate in the QED case or ii) the quarks and gluons fragment or some kind of QCD phase 
transition occurs in the QCD case. 

To simplify matters, we will assume that the plasma cools to the temperature A, at which 
point hadronization occurs. The emergent photon spectrum is a convolution of the quark 
and gluon spectrum with the pion fragmentation function 0, |5| and the Lorentz-transformed 
neutral pion decay into photons 

dN, _ f (^ +m l)/4E, dNndg^Ej 



where dg nl /dE = 2/{^m 1T f3) is the number of photons of energy E created by an pion 
decaying isotropically in its rest frame, 7 = (1 — /5 2 ) -1 , and j3 is the velocity of the pion 
(note E 7 is Doppler shifted, m n /(2 r y(l + (3)) > E 1 > 771^/(27(1 — /?))) . The pion spectrum 
is 01 

dJ% = v /■«=<» dN :i (Q,T )dg :Jn (Q,E n ) 

dE w yJ Q =E^ dQ dE v ^ U 



where dNj/dQ is the flux spectrum of the quark or gluon j of energy Q at the time of 
fragmentation, which is on the outer edge of the photosphere where the plasma is at a tem- 
perature of T = A, and dgj w /dE n is the relative number of pions with energy E n produced 
by j H- We will use dg/dE n = (15/16)2r 3 / 2 (l — z) 2 , where z = E n /Q 0. We approximate 
the quark and gluon spectrum in the observer frame by boosting a thermal spectrum at 
temperature T with the lorentz gamma factor 7 P of the outer edge of the outward moving 
photosphere 

dNj __ a iy v Q 2 f 1 (l-/3cosfl)cosfl ^ 



dQ 2n 2 Jcose=o exp (7 P o;(l — (5 cos 9)/T ) ± 1 

where aj is the number of internal degrees of freedom of particle j, the sign in the denominator 
is for fermions or bosons, and we have we integrated over the surface of the photosphere 
with radius r p . Using ref. ||, we approximate 7 P ~ (T/A) 1 / 2 and r p = 7 P /A. An accurate 
calculation of the spectrum would require using a boltzmann equation to determine the exact 
spectrum of the nearly thermal quarks and gluons, and a jet fragmentation code for the final 
decay into pions and photons. We expect the approximations used will be correct within a 
factor of order unity. 

In Figure 1 we show the photon spectrum calculated using the above formulas both 
including and excluding the QCD photosphere. Without the photosphere, the quark and 
gluon spectrum is simply a blackbody spectrum (with spin and finite size effects) at the 
temperature T. We see that the main difference between these two spectra is that the 
photosphere spectrum has many more lower energy photons. This is physically due to the 
fact that the photosphere processes high energy quarks and gluon into many lower energy 
ones, and these eventually fragment and decay into lower energy photons. We also plot the 
spectrum of photons emitted directly from the black hole which peak at approximately 5T. 
The peak height is several orders of magnitude lower than the peak from quarks and gluons, 
mostly due to the QCD degrees of freedom which can decay into photons is large |4[|. The 
total number of photons emitted from the QCD photosphere, which peaks at energies of 
about lOOMeV, is 



/ T 



N,*2x 10 24 (^— j sec" 1 (5) 

which scales as T 2 , which is a stronger function of black hole temperature than the T depen- 
dence of black body particle emission, or the T 3//2 dependence obtained in the fragmentation 
model of MacGibbon and Weber. Note that once the black hole is above the QED critical 
temperature of about 45GeV, the direct photons will also be processed through a QED photo- 
sphere and be degraded to low energies || . We have included the QED photosphere spectrum 
in figure 1, which peaks between 1 and 10 MeV and produces N y « 5x 10 28 (T/100GeV) 3 / 2 s- 1 . 
The extra power of T 1 / 2 in the QCD case comes from the multiplicity of fragmentation of 
quarks into pions. 

Let us examine several observational consequences of the QCD photosphere. The most 
important consequence involves the search for individual nearby primordial black holes. 
When developing a strategy and interpreting the results of a direct search for expiring black 
holes, one must consider the density, emergent spectrum and lifetime of the black holes, all of 
which are a function of black hole mass, and the detector sensitivity and background, which 
are a both a function of photon energy. First, one can determine the optimum energy range 



to search for these black holes by considering the background : since the observed gamma 
ray background scales approximately as E~ 2A in the range lMeV to lOGeV, one can show 
using the example of the spectrum of figure 1 that for black holes with T < lOTeV (higher 
T black holes have lifetimes shorter than Is, and emit negligible amounts of radiation), the 
optimum signal to background lies in the range of 1 to lOGeV. 

Next, one can determine the optimum black hole mass to which the detector is sensitive. 
Naturally, if Ibh is the photon emission rate of the black hole, and I is the limiting flux 
to which the detector is sensitive, then C. < Ibh/^tccI 2 , where d is the distance to the black 
hole. We know that, Ibh oc M~ 2 , however, the lifetime of the black hole r oc M 3 . One 
must, therefore, also require that the total lifetime integrated number of photons incident 
on the detector is (at least) greater than unity. To be more realistic, let us require the 
observed amount of photons iV 7 > 10. Thus, (lBH/^d 2 )Ar > 10, where A is the area of 
the detector. When both conditions are met, one finds r > 10/(£A). The optimal observing 
conditions thus occur for black holes which just meet this criterion. By using eq. |^ and 
noting that r m M 3 /3a (see e.g. 0), we can roughly estimate the maximum observing 
distance d max ~ 0.2 pc (A/2000cm 2 ) 1 / 3 (10- 9 cm- 2 S - 1 /£) 1 / 6 . 

Let us consider two satellites: EGRET J7[, which has already accumulated several years 
of data, and GLAST ||, which is still in its planning stages. In figure 2 we show the 
maximum distance each satellite can observe a small black hole, and determine that for 
EGRET (GLAST), d max ~ 0.11 (0.31) pc. 

Once one has found a limit on the maximum distance one can observe these black holes, 
one can refashion this limit into other useful limits. MacGibbon and Carr ||, and Halzen et 
al. H use the Page-Hawking limit flQ |, which is discussed below, to find a generic limit on 



the local density of black holes below a mass M to be n bh < Af((/3)(M/MA 3 pc~ 3 , where 
M* k 5x 10 14 g is the mass of a black hole which has a life of the age of the universe, 
M ~ 10 _4 pc -3 is Page-Hawking limit on the average density of M < M* black holes, and 
( is the local density enhancement of black holes compared to the universal average. This 
general scaling solution is valid up to masses M* and is related to the fact that the black hole 
lifetime is proportional to M 3 . We can make a similar limit on the value of (J\f by assuming 
that if EGRET and GLAST find no black holes after a time r = (iA)' 1 , then this will 
optimally constrain the density black holes with masses M whose lifetime is r. However, 
there is a subtle beneficial effect: after observing for a time t > r, larger mass black holes, 
which have higher number densities, will have decayed to mass Mq. Since dM/dt = —a/M 2 
(see e.g. ||), one can begin to constrain black holes with mass M 3 = 3aSt + M$, where 
St = t — t. If the constraint on the density of M black holes is no < (Qd^/S)" 1 , where Q is 
the solid angle covered by the detector, then we obtain the limit 

9 / Ml \ 1.5 x 10 10 /yr\ 

^ nd 3 \M 3 + 3aSt) ~ Qd 3 KStJ' {) 

as long as t > t. Simply put, observing for a longer time allows one to observe larger, more 
densely populated black holes which decay to the optimum observing mass M , and this 
allows for better limits on (J\f. With four years of observation and assuming no black holes 
are found with EGRET or GLAST, one can place the limits 

CAT < 5.7 x 10 12 pc~ 3 (EGRET) 



(Af < 8.5 x 10 10 pcT 3 (GLAST) (7) 

Notice that this limit scales as yi/A. As pointed out by Halzen et al. ||, luminous matter 
clusters locally by a factor on the order of 10 7 , thus since Af < 10 4 , EGRET may be able to 
see nearby primordial black holes. One can also place a limit on the local rate R at which 
black holes expire per unit volume. Roughly, R « n(M)/r(M) « (Afa/M 3 . For EGRET 
(GLAST) one finds R < 1100 (17) pcr 3 yr _1 - 

Comparison of these results to previous results is difficult because non-standard particle 
models were often used. For example, there have been several limits placed on R based on 
the Hagedorn model [p| , which presumes that the number of degrees of freedom increase 



exponentially with temperature, a fact that has no experimental support for energies up to 
about ITeV. It is important to note, however, that Alexandreas et al. [Q use the standard 



model of MacGibbon and Webber, and conclude from data of a TeV air shower array that 
R < 8.5 x 10 5 pc _3 yr _1 . However, as illustrated in figure 3, the photosphere dramatically 
alters the lifetime integrated spectrum, and the expected flux of photons above ITeV is about 
four orders of magnitude lower than the MacGibbon and Webber result, which translates 
into a limit on R six orders of magnitude weaker (higher). Note especially that when the 
photosphere is included, the lifetime integrated spectrum decreases as E~ A instead of E~ 3 . 
In fact the total number of photons above energy Ed produced by a black hole of initial 
temperature T is A" 7 « 2 x 10 u (GeV/E d ) 3 , which is valid for energies E d > (T/GeV) 1/2 GeV. 



This is to be compared to the result of Halzen et al, which finds A" 7 oc Ej . This lower flux 
at high energies renders any search for primordial black holes with TeV air shower arrays 
impractical. 

It is interesting to note that Wright |R| points out that the anisotropic component of 



the gamma ray background may be explained by the presence of primordial black holes 
clustered in the halo of our galaxy. If black holes are responsible for the anisotropy, then he 
finds (Af = (2 — 12) /h x 10 9 pc~ 3 , where h is the Hubble parameter, and this is only about 
an order of magnitude lower than the estimated detection limits of the GLAST project. 

Another consequence involves the Page-Hawking limit, which constrains (or possibly 
measures) the density of primordial black holes by comparing their expected contribution to 
the gamma ray background, for a given black hole density, to the actual observed background. 
MacGibbon and Carr || , and Halzen et al. 0, have studied this approach in detail, using a 
model of direct fragmentation of quarks and gluons (i.e. no photosphere), and they conclude 
that black holes with mass less than about 10 16 g cannot contribute more than a fraction 
flbh < 10~ 8 of the critical density of the universe. 

We have performed the same calculation including the effects of the QCD photosphere 
(see figure 3) and verified, as suggested by Heckler |§, that the QCD photosphere does not 
significantly change (i.e. less than a few percent) the limit on Vt^ found by Halzen et al. 0. 
The effect is small because the QCD photosphere becomes important only at temperatures 
above A, and black holes with these temperatures do not significantly contribute to the total 
number of photons at lOOMeV. There is one important difference in the spectrum: above 
energies of about 300MeV the spectrum including the photosphere is proportional to E~ A , 
whereas the results of Halzen et al. show a E~ 3 dependence. The steeper slope is due to 
the degrading of high energy quarks into lower energies as they are processed through the 
photosphere. 



As a final note, since the QCD photosphere will emit charged pions as well as neutral 
ones, there will also be a flux of neutrinos, electrons, and positrons up to several orders of 
magnitude larger than previously assumed, with spectra very similar to the photon spectrum 
in Figure 1. This will make constraints on (or possibilities of detection of) high energy 



neutrinos |14|, and positrons || from black holes much more important. 
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Figure 1: Instantaneous emergent spectra. Thick lines are for a M = 10 9 g (T = lOTeV) 
black hole. Thin lines for M = 10 11 g (T = lOOGeV). The solid lines are spectra which 
include photosphere: the ones peaking at about lOOMeV are the emergent spectra of QCD 
photosphere, the ones peaking at about lMeV are for the QED photosphere. For compari- 
son, we plot the dotted lines which are the direct fragmentation results of MacGibbon and 
Webber, and the dashed lines which are the direct photon emission spectra. The actual full 
spectrum is the addition of the two solid lines. 
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Figure 2: Observational distance limits. The dotted line is the minimum average distance to 
the nearest black hole of mass M, using the Page-Hawking constraint and assuming a clus- 
tering factor ( = 10 6 . The solid (dashed) lines pertain to GLAST (EGRET). Lines sloping 
upward show the maximum distance a black hole can be detected for the GLAST (EGRET) 
point source sensitivity, and the line sloping downward shows the maximum distance a black 
hole can be and still cast 10 photons on the detector, integrated over the lifetime of the hole. 
As the observation time of the detector increases, the intersection of the lines effectively 
moves to the left (see eq. ||), possibly intersecting with the Page-Hawking constraint. 
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Figure 3: Lifetime integrated emergent spectra for a M = 10 10 g (T = 1 TeV) black hole, 
which has a lifetime of about 400s. Notation is the same as in Figure 1. Notice that the 
slope of the photosphere spectra (solid line) runs as E~ 4 , whereas for the direct fragmentation 
picture dotted line) the slope runs as E~ 3 . 
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